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In the first part of the paper we show that the constraining potentials introduced to mimic 
entanglement effects in Edward's tube model and Flory's constrained junction model are diagonal in 
the generalized Rouse modes of the corresponding phantom network. As a consequence, both models 
can formally be solved exactly for arbitrary connectivity using the recently introduced constrained 
mode model. In the second part, we solve a double tube model for the confinement of long paths 
in polymer networks which is partially due to crosslinking and partially due to entanglements. Our 
model describes a non-trivial crossover between the Warner-Edwards and the Heinrich-Straube tube 
models. We present results for the macroscopic elastic properties as well as for the microscopic 
deformations including structure factors. 
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I. INTRODUCTION 

Polymer networksS are the basic structural element of 
systems as different as tire rubber and gels and have a 
wide range of technical and biological applications. From 
a macroscopic point of view, rubber-like materials have 
very distinct visco- and thermoelastic properties. □'□ They 
reversibly sustain elongations of up to 1000% with small 
strain elastic moduli which are four or five orders of mag- 
nitude smaller than for other solids. Maybe even more 
unusual are the thermoelastic properties discovered by 
Gough and Joule in the 19th century: when heated, a 
piece of rubber under a constant load contracts, and, 
conversely, heat is released during stretching. This im- 
plies that the stress induced by a deformation is mostly 
due to a decrease in entropy. The microscopic, statistical 
mechanical origin of this entropy change remained ob- 
scure until the discovery of polymeric molecules and their 
high degree of conformational flexibility in the 1930s. In 
a melt of ideutical chains polymers adopt random coil 
conformational with mean-square end-to-end distances 
proportional to their length, (r^) ~ A^. A simple sta- 
tistical mechanical argument, which only takes the con- 
nectivity of the chains into account, then suggests that 
flexible polymers react to forces on their ends as linear, 
entropic springs. The spring constant, k = ^^^y-, is pro- 
portional to the temperature. Treating a piece of rubber 
as a random network|-p£,non-interacting entropic springs 
(the phantom modela 0) qualitatively explains the ob- 
served behavior, including — to a first approximation — 
the shape of the measured stress-strain curves. 

In spite of more than sixty years of growing quali- 
tative understanding, a rigorous statistical mechanical 
treatment of polymer networks remainS|-a challenge to 
the present day. Similar to spin glasses ,lI the main dif- 
ficulty is the presence of quenched disorder over which 
thermodynamic variables rWf^ to be averaged. In the 
case of polymer networkSjETtj the vulcanization process 
leads to a simultaneous quench of two different kinds of 
disorder: (i) a random connectivity due to the introduc- 



tion of chemical crosslinks and (ii) a random topology 
due to the formation of closed loops and the mutual im- 
penetrability of the polymer backbones. Since for instan- 
taneous crosslinking monomer-monomer contacts and en- 
tanglements become quenched with a probability propor- 
tional to their occurrence in the melt, ensemble averages 
of static expectation values for the chain structure etc. 
are not affected by the vulcanization as long as the sys- 
tem remains in its state of preparation. 

For a given connectivity the phantom model Hamilto- 
nian for non-interactins, polymer chains formally takes a 
simple quadratic form,m3 so that one can at least for- 
mulate theories which take the rjaHdam connectivity of 
the networks fully into account .llilil^l The situation is 
less clear for entanglements or topological constraints, 
since they do not enter the Hamiltonian as such, but 
divide phase space into accessible and inaccessible re- 
gions. In simple cases, entanglements can be character- 
ized bT|t,[i.opological invariants from mathematical knot 
theory.Ba However, attempts to formulate topeiogical 
theories of rubber elasticity (for references seeEj) en- 
counter serious difficulties. Most theories therefore omit 
such a detailed description in favor of a mean-field ansatz 
where the different parts of the network are thought to 
move in a deformation-dependent elastic matrix which 
exerts restoring forces towards some rest positions. These 
restoring forces may be due to chemical crosslinks which, 
localize random paths through the network in spaceliEl 
or to entaaglernents. The classical theories of rubber 
elasticityEnlTo assume that entanglements act only on 
the. |Crcisslinks or junction points, while the tube mod- 
clsH'El tj stress the importance of the topological con- 
straints acting along the contour of strands exceeding a 
minimum "entanglement length" , A^g. Originally devised 
for polymer networks, the tube concept is particularly 
successful in explaining the extremely long relaxation 
times in non-crosslinked polymer melts as the result ofj-a 
one-dimensional, curvilinear diffusion called reptationEII 
of linear chains of length N ^ within and finally out 
of their original -tubes. Over the laat_decade computer 
simulationgLjc3~E3 and experimentsElTES have finally also 
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collected mounting evidence for the importance and cor- 
rectness of the tube concept in the description of polymer 
networks. 

More than thirty years after its introduction and 
in spite of its intuitivity and its success in provid- 
ing a pUBi fie j d .view on entangled polymer networks and 
mcltSjHS'EaE^ there exists to date no complete solution 
of the Edwards tube model for polymer networks. Some 
of the open problems are apparent from a ijeceot con- 
troversy on the interpretation of SANS data.E3~0 Such 
data constitute an important experimental test of the 
tube concept, since they contain information on the de- 
gree and deformation dependence of the confinement of 
the microscopic chain motion and therefore allow for a 
more detailed test pf theories of rubber elasticity than 
rheological data.E3E3 

On the theoreticaL-jSide, the original approach of 
Warner and Edwardstljised mathematically rather in- 
volved replica methodscj to describe the localization of 
a long polymer chain in space due to crosslinking. The 
replica method allows for a very elegant, self-consistent 
introduction of constraining potentials, which confine in- 
dividual polymer strands to random- walk like tubular re- 
gions in space while ensemble averages over all polymers 
remain identical to thpsp of unconstrained chains. Later 
Heinrich and StraubeE3c3 recalculated these results for 
a solely entangled system where they argued that there 
are qualitative differences between confinement due to 
entanglements and confinement due to crosslinking. In 
particular, they argued that the strength of the confin- 
ing potential should vary afhnely with the macroscopic 
strain, resulting in fluctuations perpendicular to the tube 
axis which vary only like the square root of the macro- 
scopic strain. 

Replica calculations provide limited insight into phys- 
ical mechanism aad make approximations which are dif- 
ficult to control.E3 It is therefore interesting to note 
that Flory was able to solve the_in many respects sim- 
ilar, constrained-junction modeltO without using such 
methods. Recent refinements of the constrainpi-junction 
model such as the constrained-chain modeH and the 
diffused-constrained modelEa have more or less converged 
to the (Heinrich and Straube) tube model, even though 
the term is not mentioned explicitly. Another variant of 
this mOjdel was recently solved by Rubinstein and Pa- 
nyukov.CJ In particular, the authors illustrated how non- 
trivial, sub-affine deformations of the polymer strands 
result from an affinely deforming confining potential. 

While tube models are usually formulated and dis- 
cussed in real space, two other recent papers have pointed 
independently to considerable simplifications_pf the cal- 
culations in mode space. Read and McLeishEEl were able 
to rederive the Warner-Edwards result in a particularly 
simple and transparent manner by showing that a har- 
monic tube potential is diagonal in the Rouse modes of a 
linear chain. Complementary, one of the present authcps 
introduced a general constrained mode model (CMM) ,63 
where confinement is modeled by deformation dependent 



linear forces coupled to (approximate) eigenmodes of the 
phantom network instead of a tube-like potential in real 
space. This model can easily be solved exactly and is 
particularly suited for the analysis of simulation data, 
where its parameters, the degrees of confinement for all 
considered modes, are directly measurable. Simulations 
of defect-free model polymer networks^ under strain ana- 
lyzed in the framework of the CMMIiJ provide evidence 
that it is indeed possible to predict macroscopic restor- 
ing forces and microscopic deformations from constrained 
fluctuation theories- In particular, the rssults support 
the choice of Flory,lllLHeinrich and Straube,ca and Rubin- 
stein and PanyukovEJ for the deformation dependence of 
the confining potential. In spite of this success, the CMM 
in its original form suffers from two important deficits: 

(i) due to the multitude of independent parameters it is 
completely useless for a comparison to experiment and 

(ii) apart |iom recovering the tube model on a scaling 
level, Refo remained fairly vague on the exact relation 
between the approximations made by the Edwards tube 
model and the CMM respectively. 

In the present paper we show that the two models are, 
in fact, equivalent. The proof, presented in section [IB 



is a generalization of the result by Read and McLeish 
to arbitrary connectivity. It prpijides the liiik between, 
the considerations of EichingerjHl Graessley,c3 Mark,c3 
and others on the dynamics of (micro) phantom networks 
and the ideas of Edwards and Flory on the suppression 
of fluctuations due to entanglements. As a consequence, 
the CMM can be used to formally solve the Edwards tube 
model exactly, while in turn the independent parameters 
of the CMM are obtained as a function of a single param- 
eter: the strength of the tube potential. Quite interest- 
ingly, it turns out that the entanglement contribution to 
the shear modulus depends on the connectivity of the net- 
work. In order to explore the consequences, we discuss in 
the second part the introduction of entanglement effects 
into the Warner-Edwards model, which represents the 
network as an ensemble of independent long paths com- 
prising many strands. Besides recovering some results by 
Rubinstein and Panyukov for entanglement dominated 



Je chain structure fac- 
Finally we propose a 



systems, we also calculate the^ 
tor for this controversial case.r^'~ 
"double tube" model to describe systems where the con- 
finement of the fluctuations due to crosslinks and due to 
entanglements is of similar importance and where both 
effects are treated within the same formalism. 



II. CONSTRAINED FLUCTUATIONS IN 
NETWORKS OF ARBITRARY CONNECTIVITY 

A. The Phantom Model 

The Hamiltonian of the phantom model0~l is given by 
"^ph = ^J2'{ij) ^ij^ where {i,j) denotes a pair of nodes 
i,j € 1, . . . , M which are connected by a polymer chain 
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acting as an entropic spring of strength k 



3k bT 



and 

, (t) — ri (t) — Tj (t) the distance between them. In or- 
der to simphfy the notation, we always assume that all 
elementary springs have the same strength k. The prob- 
lem is most conveniently studied using periodic boundaot 
conditions, which span the network over a fixed volumetJ 
and define the equilibrium position Ri = {Xi,Yi, Zi). A 
conformation of a network of harmonic springs can be 
analyzed in terms of either the bead positions fi{t) or 
the deviations Ui (t) of the nodes from their equilibrium 
positions Ri. In this representation, the Hamiltonian sep- 
arates into two independent contributions from the equi- 
librium extensions of the springs and the fluctuations. 
For the following consideratioas it is useful to write fluc- 
tuations as a quadratic form.EJ Finally, we note that the 
problem separates in Cartesian coordinates a = x,y, z 
due to the linearity of the springs. In the following we 
simplify the notation by writing the equations only for 
one spatial dimension: 
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Here u denotes a M-dimensional vector with {u)i = 
(ui)x- K_ is the connectivity or Kirchhoff matrix whose 
diagonal elements {K_)ii = fik are given by the node's 
functionality (e.g. a node which is part of a linear chain 
is connected to its two neighbors so that fi — 2 in con- 
trast to a four- functional crosslink with fi — 4). The 
off-diagonal elements of the Kirchhoff matrix are given 
by {K_)ij — —k, if nodes i and j are connected and by 
(K )ij = otherwise. Furthermore, we have assumed 
that all network strands have the same length. 

The fluctuations can be written as a sum over indepen- 
dent modes Cp which are the eigenvectors of the Kirch- 
hoff matrix: K_ep = kpCp where the Cp can be chosen 



to be orthonormal e p • e p/ 



i5pp/. The transformation to 



the eigenvector representation u = S_u and back to the 
node representation u = S_~^u is mediated by a matrix 
S_ whose column vectors correspond to the Cp. By con- 
struction, 5 is orthogonal with 5 * —S_^^. Furthermore, 
the Kirchhoff matrix is diagonal in the eigenvector repre- 



sentation [K_ )pp 
then reduces to 



{S_ ^K_S_)pp = kp. The Hamiltonian 
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Since the connectivity is the result of a random process, it 
is difficult to discuss the properties of the Kiiphhoff ma- 
trix and the eigenmode s^ctrum in general.ElrCj The fol- 
lowing simple arguments ignores these difficulties. The 
idea is to relate the mean square equilibrium distances 
{Xfj) to the thermal fluctuations of the phantom net- 
work. 

Consider the network strands before and after the 
formation of the network by end-linking. In the melt 



state, the typical mean square extension ^f^) is en- 
tirely due to thermal fluctuations, while (r) = 0. In 
the crosslinked state, the strands show reduced ther- 
mal fluctuations (u^^) around quenched, non- vanishing 

mean extensions {Rfj)- However, the ensemble aver- 

- (ufA is not affected 



age of the total extension (Rij) + { 
by the end-linking procedure. The ffuctuation contribu- 



tion 



{ufj) depends on the connectivity of the network 



and can be estimated using the equipartition theorem. 
The total thermal energy in the fluctuations, Ufiuc, is 
given by ^ksT times the number of modes and therefore 



}kBTN„a 



J ^ksTN strands, whcrC Nnodes 



Ufluc = ^KB 

and N strands are the number of junction points and net 
work strands, which are related by N strands — 2-^nodes 
in an /-functional network. EquatiH|g-|the thermal energy 



per mode to ^{uf,), one obtainsoc 
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(3) 
(4) 



Using these results, one can finally estimate the elas- 
tic properties of randomly cross- or endlinked phantom 
networks. Since the fluctuations are independent of size 
and shape of the network, they do not contribute to 
the elastic response. The equilibrium positions of the 
junction points, on the other hand, change affinely in 
the macroscopic strain. The elastic free energy density 
due to a volume-conserving, uni-axial elongation with 



All = A 



-1/2 



A is simply given by: 



AFp^(A) = A^ + 



A 



2 ) 

strand' 



Pstrand 



1 ^ y ) Pstrand , (5) 



where Pstrand is the number density of elastically active 
strands. For incompressible materials such as rubber, the 
shear modulus is given by ^ of the second derivative of 
the corresponding free energy density with respect to the 
strain parameter A. In response to a finite strain, the 
system develops a normal tension <tt'- 



Gph 



1 1 <fAFph{\) 
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A=l 



= I 1 — J ) PstrandksT 



OT = (A^ - -^G, 



(6) 

(7) 

(8) 



Experimentally observed stress-strain curves show devia- 
tions from eq ^. Usually the results are normalized to the 
classical prediction and plotted versus the inverse strain 
1/A, since they often follow the semi-empirical Mooney- 
Rivlin form 
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(9) 



C. Solution and Disorder Averages: The 
Constrained Mode Model (CMM) 



Since the total Hamiltonian of the CMM 



B. The Constraint Hamiltonian 

Most theories introduce the entanglement effects as 
additional, single-node terms into the phantom model 
Hamiltonian, which constrain the movement of the 
monomers and junction points. The standard choice are 
anisotropic, harmonic springs of strength / (A) between 
the nodes and points CiW which are fixed in space: 



-He 



(10) 



While all models assume that the tube position changes 
affinely with the macroscopic deformation. 



e,(A) = AC. (A = 1) 



(11) 



there are two different choices for the deformation depen- 
dence of the confining potential: 



/ (A) = ? (A = 1) 
T(A) = a"^T(A = 1). 



(12) 
(13) 



Since this choice of Ticonstr leaves the different spatial 
dimensions uncoupled, we consider the problem again in 
one dimension and express Ticonstr in the eigenvector rep- 
resentation of the Kirchhoff matrix of thejjnconstrained 
network. Using v{X) — (A) — ^(A) = Aw (A — 1) one 
obtains 

T-Lconstr = "Y^di - nY {^L " H) 

= M(^-i^ -g-'vYii-'u-i-'v) 

= ^-^in~v)%g-Hu~v) (14) 
= ^(M -£)*(£-£) 

E?(A) _ 2 

p 

Thus the introduction of the single node springs does not 
change the eigenueciors of the original Kirchhoff matrix. 
The derivation of eq |lj, which is the Hamiltonian of the 
Constrained Mode Model (CMM),o is a central result of 
this work. It xurovides the liiik between the considerations 
of EichingerjHl Graessley,L£l Mark,c£l and others on the 
dynamics of (micro) phantom networks and the ideas of 
Edwards and Flory on the suppression of fluctuations due 
to entanglements. 
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(15) 



is diagonal and quadratic in the modes, both the exact 
solution of the model for given Vp and the subsequent 
calculation of averages over the quenched Gaussian dis- 
order in the Vp are extremely simple.O In the following 
we summarize the results and give general expressions 
for quantities of physical interest such as shear moduli, 
stress-strain relations, and microscopic deformations. 

Consider an arbitrary mode Up of the polymer net- 
work. Under the influence of the constraining potential, 
each Cartesian component a will fluctuate around a non- 
vanishing mean excitation Up with 



C/pa(A) 



;a(A) 



kp -l- ?a(A) 



.(A). 



(16) 



Using the notation 5up = Up — Up, the Hamiltonian for 
this mode reads 



HpM ^^U'{\) + ^ (C/p„(A) - Vp^iX)) 



2 

laiX) + kp 



(17) 



Expectation values are calculated by averaging over both 
the thermal and the static fluctuations, which are due to 
the quenched topological disorder:^ 



(Ap(A)) = J dvp J dSup Ap[vp,6up] P{vpa)P{5upa). 

(18) 

Both distributions are Gaussian and their widths 

keT 



^pal 



kp -f la (A) 
kp ^ y 



(19) 
(20) 



follow from the Hamiltonian and the condition that the 
random introduction of topological constraints on the dy- 
namics does not affect static expectation values in the 
state of preparation. In particular. 



("pj = (^O + {Up 



2 > KiS 



knT 



(21) 



Eq ^ relates the strength I of the confining potential 
to the width of P{vpa). The result, (w^^) = -^^^ 



^In order to simplify the notation, we use I = /(A 
(up„) = (^^^(A = 1)) etc. 



1), 
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(f^po) = 7p^, {SuU = (1 - 7p)^ can be expressed 
conveniently using a parameter 



0<7p 



I 



< 1, 



(22) 



which measures the degree of confinement of the modes. 
As a result, one obtains for the mean square static exci- 
tations: 



kp + la ('^) / I 



kr) 



(23) 



Quantities of physical interest are typically sums over 
the eigenmodes of the Kirchhoff matrix. For example, 
the tube diameter is defined as the average width of the 
thermal fluctuations of the nodes: 



In particular, 



4, 



ksT 
Ml 



E 

p 



(24) 



(25) 



More generally, distances between any two monomers 
Tnma — fna ~ ''met in real spacc are given by 



(26) 



For the discussion of the elastic properties of the dif- 
ferent tube models it turns out to be useful to define the 



sum 




Using 
tensio: 

(7T 



- 1 j • (27) 

the confinement contribution to the normal 
Hand the shear modulus can be written as 



= (A2-l)g(A) + (l-A-i).g(A-i/2) (28) 

Gcor^str = .g(l). (29) 



D. Model A: Deformation independent strength! of 
the Confining Potential 

In order to completely define the model, one needs to 
specify the deformation dependence of the confining po- 
tential. One plausible choice is 



;^(A)= Ia(\ = \) , 



(30) 



i.e. a confining potential whose strength is strain inde- 
pendent. The following discussion will make clear, that 
this choice leads to a situation which mathematically re- 
sembles the phantom model without constraints. 

Using eq ^ the thermal fluctuations (and therefore 
also the tube diameter eq ^) are deformation indepen- 
dent and remain isotropic in strained systems. The mean 
excitations, on the other hand, vary affinely with the 
macroscopic strain. This leads to the following relation 
for the deformation dependence of the total excitation of 
the modes: 



pa 



>(A) 



1 + (Ai - 1) 



Ia 



Using eqs 



(■Wpa) rup 

I to |29| one obtains via 



kn + 1 A 



9a{\) 



E 



V ^\kp + lA 



a classical stress-strain relation: 

ariX) = {X^-X-')Ga 
keT 



Ga 



V 



E>- 

p 



(31) 
(32) 

(33) 
(34) 



E. Model B: Affine deformation of the confining 
potential 



The ansatz 



^s(A) = A Zb(A = 1) 



(35) 



goes back to Ronca and AUegraEEl and was used by Flooi. 
Heinrich and StraubeEj and Rubinstein and Panyukov.L3 
It corresponds to afflnely deforming cavities and leads 
to a more complex behavior including corrections to the 
classically predicted stress-strain behavior. 

Using eq 35, the mean excitations of partially frozen 
modes as well as the thermal fluctuations, become de- 
formation dependent. The total excitation of a mode is 
given by: 



<)(A) 



1 + (A^ - 1) 



pet/ 



kp + IsiX) 



(36) 



Only in the limit of completely frozen modes, 7^ ^ 1, one 
finds affine deformations with Upa{X) = XaUpa{X = 1). 

takes the form 



Concerning the elastic properties, eq |27 



9b{X) 



keT 
V 



E 

p 



Ib{X) 



kp + 1 six) 



while the the shear modulus can be written as 
keT 



G 



B 



V 



■E7|- 



(37) 



(38) 
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Note the different functional form of eqs |3J and 38 , Since 
< 7 < 1, the contribution of confined modes to the 
elastic response is stronger in Model A than in Model 
B. Furthermore, within Model B the interplay between 
the network connectivity (represented by the eigenmode 
spectrum {fcp} of the Kirchhoff matrix) and the confining 
potential I is different for the shear modulus eq ^ and 
the tube diameter eq E5l 



F. Model C: Simultaneous presence of both types of 
confinement 

Finally, we can discuss a situation where confinement 
effects of type A and B are present simultaneously. Cou- 
pling each node to two extra springs /^(A) = I a and 
^_Ba(A) — Isa/^a leads to the following Hamiltonian in 
the eigenmode representation: 



Ia 



{Up - VApiX)Y 



IbW 



{up - WBp(A))^. 

(39) 



Model A and Model B are recovered by setting I a re- 
spectively Ib equal to zero. Furthermore we assume, that 
both types of confinement can be activated and deacti- 
vated independently. This requires, 



LAKp 
, 2 X _ + kp 

v^Bpl ~ — ; — ■ 

iBKp 



(40) 
(41) 



In the presence of both types of confinement, the mean 
excitation of the modes is given by 



(C/p>(A) = 



{lA + lB{X)+kp)^ 

21aIb{X) {vApVBp){X) 

{lA + lB{X) + kpf 

while the thermal fluctuations are reduced to 



(42) 



(43) 



Ia + IbW + kp ' 

Finally the condition that the simultaneous presence of 
both constraints does not affect ensemble averages in the 
state of preparation requires 

/ \ ksT 

{vApVBp} = —r— ■ (44) 

Kp 

From eqs one can calculate the deformation de- 

pendent total excitation of the modes: 



1 + {X'-1) 
{kp + lA+^i 



(45) 



so that 



f\\ ksT ^ 



kp + lA+'^ {kp + lA+'-^f 



(46) 



In the present case, the shear modulus can be written as 
kBT 7Ap(l - iBp) 



Gc 



V ^ l--iAplBp 



7Sp(l - lAp) iciBp(X - lAp) + 7Ap(l - 7Sp)) 



(1 - lAplBpY 



(47) 



Note that the shear modulus is not simply the sum of 
the contributions from the A and B confinements. While 



eqs |3J and 38 are reproduced in the limits ^Bp — and 
"fAp = respectively, eq ^ reflects the fact that a mode 
can never contribute more than ksT to the shear mod- 
ulus. Thus, for 7Bp = 1 (respectively 7Ap = 1) the pth 
mode contributes this maximum amount independently 
of the value of ^Ap (respectively ^Bp)- 

An important point, which holds for all three models, 
is that it is not possible to estimate the confinement con- 
tribution to the shear modulus from the knowledge of the 
absolute strength Ia,Ib of the confining potentials alone. 
Required is rather the knowledge of the relative strengths 
lAp, iBp which in turn are functions of the network con- 
nectivity. 



G. Discussion 

It is not a priori clear, whether entanglement effects 
are more appropriately described by model A or model 
B. While model A has the beaefit of simplicity, Ronca and 
Allegra proposed model B,Ea because it leads (on length 
scales beyond the tube diameter) to the conservation of 
intermolecular contacts under strain. Similar conclusions 
were drawn bypHeinrich and Straubeo and Rubinstein 
and Panyukov.EJ In the end, this problem will have to 
be resolved by a derivation of the tube model from more 
fundamental topological considerations. For the time be- 
ing, an empirical approach seems to be the safest option. 
Fortunately, tli©,evidence provided by experimentsE3 and 
by simulationalJ points into the same direction. 

Since details of the interpretation of the relevan t ex- 
periments are still controversial (see Section III D 3 ) , we 
concentrate on simulation results where the strain depen- 
dence of approximate eigenmodes of the phantom model 
was measured directly.o Figure shows a comparison of 
data obtained for defect-free model polymer networks to 
the predictions eq ^ of Model A and eq ^ of Model 
B. The result is unanimous. We therefore believe eq |35| 
and Model B to be the appropriate choice for modeling 
confinement due to entanglements. The shear npdulus 
of an entangled network should thus be given hyx-i 



6 



G = G 



ph 



V 



(48) 



the various 7p are no 



where in contrast to ReferenC' 

longer free parameters but depend through eqM on a sin- 
gle parameter: the strength I of the confining potential, 
which is assumed to be homogeneous for all monomers. 
The difficulty of this formal solution of the generalized 
constrained fluctuation model for polymer networks is 
hidden in the use of the generalized Rouse modes of the 
phantom model, wliicb. are difficult to obtain for real- 
istic connectivities 00 A useful ansatz for end-linked 
networks is a separation into independent Flory-Einstein 
respectivelji,Rouse modes for the crosslinks and network 
strands. Ejo In fact, the simulation results presented in 
Figure |^ are based on such a decomposition. 

For randomly crosslinked networks with a typically 
exponential strand length polydispersity, the separation 
into Flory-Einstein and single-chain Rouse modes ceases 
to be useful. In this case, we can think of two radically 
different strategies: 

• To keep the network connectivity in the analysis. 
For example, there is no principle reasan, why the 
methods presented by Sommer et al.Ell and Ever- 
aeralil could not be combined, in order to investi- 
gate the strain dependence of constrained general- 
ized Rouse modes in computer simulations. Note, 
however, that this completely destroys the self- 
averagingjiroperties of the approximation used in 
ReferenceEJ. Analytic progress in the evaluation of, 
for example, eq 38 for the entanglement contribu- 



tion to the shear modulus requires information on 
the statistical properties of the eigenvalue spectra 
of networks generated by random crosslinking. To 
our knowledge, the only available results were ob- 
tained numerically by Shy and Eichinger.113 Note, 
that model C is irrelevant, if one is able to carry out 
calculations with the proper network eigenmodes. 

• To average out the connectiwity effects in tube mod- 
els for polymer networks. t3 In the second part of 
the paper, we will consider linear chains under the 
influence of two types of conflnement: network con- 
nectivity and entanglements. 



III. TUBE MODELS 

In SANS experiments of dense polymer melts, it is pos- 
sible to measure sirLgle chain properties by deuterating 
part of the polymers.c3 If such a system is flrst crosslinked 
into a network and subsequently subjected to a macro- 
scopic strain, one can obtain information on the micro- 
scopic deformations of labeled random paths through the 
network.CJ In order to interpret the results, they need 



to be compared to the predictions of theories of rubber 
elasticity. Unfortunately, for randomly crosslinked net- 
works it is quite difficult to calculate|±he,pJevant struc- 
ture factors even in the simplest cases.ESE^O Because the 
crosslink positions on different precuESor chains should 
be uncorrelated, Warner and EdwardsE3 had the idea to 
consider a tube model, where the crosslinking effect is 
"smeared out" along the chain. To model confinement 
due to crosslinking, they used (in our notation) Model 
A, since this ansatz reproduces the essential properties of 
phantom models (affine deformation of equilibrium posi- 
tions and deformation independeacc of fluctuations). In 
contrast, Jieinrich and StraubecJ, and Rubinstein and 
PanyukovH treated conflnement due to entanglements 
using Model B. Obviously, both effects are present si- 
multaneously in polymer networks. In the following, we 
will develop the idea that in order to preserve the qualita- 
tively different deformation dependence of the two types 
of confinement, they should be treated in a "double tube" 
model based on our Model C. 

Before entering into a detailed discussion, we would 
like to point out a possible source of confusion related to 
the ambiguous use of the term "tube" in the literature 
(including the present paper). A real tube is a hollow, 
cylindrical object, suggesting that in the present context 
the term should be reserved for the confining potential 
described by quantities such as ^i,Vp, I. It is in this sense 
that we speak of an "affinely deforming tube" . However, 
a harmonic confining tube potential is a theoretical con- 
struction which is difficult to visualize. For example, in 
the continuum chain limit used below, the forces exerted 
"per monomer" become infinitely small corresponding to 

— > 00, ^ 0. On the other hand, the term tube is 
often associated with the tube "contents" , i.e. the super- 
position of the accessible polymer configurations charac- 
terized via a locally smooth tube axis (the equilibrium 
positions Up) and a tube diameter dr (defined via the 
fiuctuations Sup). JThis second definition refers to mea- 
surable quantitiesJlj Which kind of tube we are referring 
to, will hopefully always be clear from the context and 
the mathematical definition of the objects under discus- 
sion. 

In the case of linear polymers, the phantom model re- 
duces to the Rouse model with vanishing equilibrium po- 
sitions Ri = 0. As a consequence, there are no strain 
effects other than those caused by the confinement of 
thermal fluctuations. In particular, the "intrinsic" phan- 
tom modulus vanishes (see eq |]). Since the networks are 
modeled as superpositions of independent linear paths, 
we have to introduce conflnement of type A in order to 
recover the phantom network shear modulus Gph in the 
absence of entanglements. 

In the Rouse model the Kirchhoff matrix takes the sim- 
ple tridiagonal form 
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A. The Warner-Edwards Model 



K = k 



/-2 
1 



Vo ••• 



1 -2/ 



(49) 



and, depending on the boundary conditions, is diago- 
nalized by transforming to sin or cos modes using the 
transformation matrix 



OP 



1 /. ip 

— exp ITT — 

N \ N 



(50) 



The eigenvalues of the diagonahzed Kirchhoff matrix 
{K)pp = {S~^KS)pp = kp are given by 



fc„ = 4fcsin2 (P^) 
^ \2nJ 



(51) 



If we consider a path with given radius of gyration i?^, 
the basic spring constant is given by fc = ^2R?^ • ^^'^ 
continuous chain limit {N oo), sums over eigenmodes 
can be approximated by integrals. For example, one ob- 
tains from eq an expression for the tube diameter 



Warner and Edwardslij used the replica method to cal- 
culate the conformational statistics of long paths through 
randomly crosslinked phantom networks. The basic idea 
was to represent the localization of the paths in space due 
to their integration into a network by a coarse-grained 
tube-like potential. Recently, it was shown by Read 
and McLeishoE3 that the same result could be obtained 
along the lines of the following, much simpler calculation, 
where we evaluate Model A for linear polymers. 

Evaluation of the integrals eqs ^ and ^ yields for the 
deformation independent tube diameter and the internal 
distances: 



^2 

"■TAc 



'a(A)) 



2RI 



2y/kU 
A^lx — 



(57) 
(58) 



+ (1- A 



2\ "-TAa 



1 - e 



We note that the latter equation can be rewritten in the 
form 



4a 



1 r 

NJo 

keT 



dp 



hp -\- 1 



keT 



y/l {ik + I) 2Vlk 



(52) 



which could be further simplified, since in this limit the 
springs representing a chain segment between two nodes 
are much stronger than the springs realizing the tube, i.e. 
fc > /. 

For normally distributed internal distances fxx' be- 
tween points X — jj^, x' — on the chain contour the 
structure factor is given by 



S{q,X)= / dx dx' exp 

Jx=0 Jx'=0 



1 ^ 

2 ^ 

Q = l 



'Z^i^L'JA)) 



(53) 



In the present case, eq Eq reduces to 



2 

XX' OL 



dv{ulj\)) 



^lirpx 



^lirpx' 



In the undeformed state, 

(rL-JA = l)) =2i?>-x'| 



(54) 



(55) 



so that the structure factor is given by the Debye func- 
tion: 



^(<f,A = l) 



2N 



expi~q^Rl)-l + q^Rl) . (56) 



.'o(A)) 



'Jl)> 



(A^-1)(^L'J1)> 



Rl\x^ 



d"^ 

"■TA 



(59) 



fAiy) = 1 + 



y 



with a universal scaling function fA(jj) which does not de- 
pend explicitly on the deformation, eq |5^ measures the 
degree of afHneness of deformations on different length 
scales. Locally, i.e. for distances inside the tube with 
Rg\x — x'\ ^ corresponding to y ^ 1, the polymer 
remains undeformed. Thus limy_>o/(j/) = 0. Deforma- 
tions become affine ioi R^\x — x'\ ^ d"^ and y ^ 1, where 
f{y) tends to one. 

Furthermore, one obtains for the shear modulus and 
the stress-strain relation 



9a{X) 
Ga 



pb'^^/klA 
6 

1 plPknT 

4 4a 
1, 



^Ga 



ariX) - (A^ - -)Ga- 



(60) 
(61) 
(62) 



so that the Mooney-Rivlin parameters are simply given 
by 



2Ci = Ga 
2C2 = . 



(63) 
(64) 
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B. The Heinrich-Straube / 
Rubinstein— Panyukov-Model 

Heinpkh and Straube,@ and Rubinstein and Pa- 
nyukovc3 have carried out analogous considerations for 
Model B, i.e. an affinely deforming tube. The relation 
between the strength of the springs I b and the tube di- 
ameter in the unstrained state is identical to the previous 
case. However, the tube diameter now becomes deforma- 
tion dependent: 



(65) 



2d = Gph + 'jGs 
2C2 = ■ 



(71) 
(72) 



Note that eq ^ holds only for A « 1. For large 
compression or extension the approximation k 3> 1{X) 
breaks djOwn and one regains the result of Heinrich and 
Straube:Ea 



or (A) 



A - ) Gb (A< 1, A> 1) . (73) 



Thus the typical width of the fluctuations changes only 
with the square root of the width of the confining poten- 
tial. Using equations eqs and one obtains for the 
mean square internal distances: 



(^L'a(A)) _ ^2 I /| 

2RI -^"l^-^l 



+ 2 



(66) 



Again, we can rewrite this result in terms of a universal 
scaling function for the degree of afhneness of the polymer 
deformation: 



(A^-1)(^L'J1)) 




1 3 p^y — 1 

2 ^ y 



(67) 



Eq ^ shows that Straube's conjectureE3~ii /a(2/) = 
fB{y) is incorrect. However, the two functions are qual- 
itatively very similar. 

For the shear modulus and the stress-strain relation we 
find 



Ipb^k 



o-t(A) 



1 ph 



TB 



Gb 

T 



j)Gb 



(68) 



(69) 



(70) 



in agreement with Rubinstein and Panyukov. In order 
to account for the network contribution to the shear mod- 
ulus, these authors add the phantom network results to 
eqs |6^, This leads to the following relations for the 
Mooney-Rivlin parameters:E3 



C. The "double tube" model 

In the following we discuss a combination of two dif- 
ferent constraints, one representing the network (Model 
A) and therefore deformation independent and the other 
representing the entanglements (Model B). Thus we use 
Model C to combine the Warner- Edwards model with the 
Heinrich-Straube /Rubinstein-Panyukov model. 

Evaluating eq Ga one obtains for the tube diameter 



1 



1 



1 



4ca(A) 



rf4 



4i3o(A) 



(74) 



The deformation dependent internal distances are 
given by: 



2RI 



Ailx- 



^'l^c^A)^- 



'4sJA)' 



— 



(A^ 



4c JA) 



i?2 



(A) {d 

TAo 



f 4sa(A)) 



4bq(A) 4ac 



(75) 



In this case, it is not possible to rewrite the result in terms 
of a universal scaling function, because the relative im- 
portance of the two types of confinement is deformation 
dependent. Introducing $( A) = 4ca(A)/4Ba(A)i eq|75| 
can be rewritten as 



(A^-1)(''L'.(1)> 
U (y) + ^X) Ub (y) 



■,$(A) 



(76) 



For the elastic properties of the double tube model we 
find: 



r, (W ^g|(A) -ryA 

y/4g%{X)+gl 
ariX) = (A^ - 1) gciX) + (l - A-^) gciX-'^') • (78) 
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Again, eq 78 only holds for moderate strains. Shear mod- 



ulus and the Mooney-Rivlin parameters are given by: 



Gc 



2Ci = 



2G| 



&G%G\ + 4G; 



(4G| 



G\ 



,3/2 



2Co = 



AG% 



(4G| 



G\f" 



(79) 
(80) 
(81) 



D. Comparison of the different tube models 

In the following we compare the predictions of the dif- 
ferent models for the microscopic deformations and the 
macroscopic elastic properties from two different points 
of view: 

1. As a function of the network connectivity, i.e. the 
ratio of the average strand length Nc between 
crosslinks to the melt entanglement length N^- For 
this purpose, we identify Ga with the shear modu- 
lus of the corresponding phantom network Gph'- 



1 ph'^ksT 
4 4a 



— Ga — Gph — 



2 



Nr. 



'2 - ^b^Nr 



rl^ — 
"'TA — 



4/ 



(82) 
(83) 



where we use / = 4 for our plots. Similarly, we 
choose for Gb a value of the order of the melt 
plateau modulus Ge- 



1. Elastic properties 

Figure ^ shows the shear modulus dependence on the 
ratio of the network strand length Nc to the melt entan- 
glement length Nf.. As expected Gc crosses over from 
Gph for short strands to Ge in the limit of infinite strand 
length. For comparison we have also included the predic- 
tion of Rubinstein and Panyukov, Gph + Ge- The shear 
moduli predicted by our ansatz are always smaller than 
this sum. In particular, we find G — Gph for N^ <C A^e- 
The physical reason is that in a highly crosslinked net- 
work the typical fluctuations are much smaller than the 
melt tube diameter. As a consequence, the network does 
not feel the additional confinement and the entangle- 
ments do not contribute to the elastic response. Figure ^ 
shows analogous results for the Mooney-Rivlin Parame- 
ters Gi and G2 again in comparison to the predictions of 
Rubinstein and Panyukov. Note, that G2 is not predicted 
to be strand length independent. 

Figure ^ shows the reduced force in the Mooney-Rivlin 
representation for different entanglement contributions 
$ to the confinement. For moderate elongations up to 
A w 2 the curves are well represented by the Mooney- 
Rivlin form. For a given shear modulus, Gi and G2 are a 
function of the entanglement contribution <I> to the con- 
finement: 



$2 



2Gi = Gc 1 



4-2$ 



2G2 = Gc 



$2 



4-2$ 



(89) 



2. The tube diameter 



1 ph^keT 

8 ^TB 



Gi 



"-TB 



Ge 
1. 

6 



3pb^kBT 

4 Ne 



-h^Ne 



(84) 
(85) 



Since eq ^ can be written in the form 
1 pPkBT 



Gc 



(2-$) 



(90) 



Assuming that the system is characterized by a cer- 
tain tube diameter (Itc or shear modulus Gc, we 
discuss its response to a deformation as a func- 
tion of the relative importance < $ < 1 of the 
crosslink and the entanglement contribution to the 
confinement: 



1-0,= ^ 



where $ is of the order 1 



(86) 
(87) 



a plot of dj^^ versus Ne/Ne looks very similar to Figure]^. 

The deformation dependence of the tube diameter 
(Figure ||) takes the form: 



lim 

A — ^oc 



"■TCa 

4cJA) 



A— *oo dxci. 



(1 - $ 

-1/4 



1/4 



(91) 
(92) 
(93) 



In the parallel direction, the entanglement contribu- 
tion to the confinement vanishes for large A so that 
lim;y_,oo d-TcwW — dTA\\- On the other hand, the entan- 
glements become relatively stronger in the perpendicular 
direction with liniA^oo dTc±W = c'ts±(A). 
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3. Microscopic deformations and structure functions 

Figure ^ compares the universal scaling functions of 
the Warner-Edwards and Heinrich-Straube/Rubinstein- 
Panyukov model defined by eqs |9[ 

More important for the actual microscopic deforma- 
tions than the difference between these two functions is 
the fact, that the distances are scaled with the deforma- 
tion dependent tube diameter. As a consequence, defor- 
mations parallel to the elongation are smaller in Model 
B than in Model A, while the situation is reversed in 
the perpendicular direction. In the general case (eq 
of Model C), the results are further complicated by the 
deformation dependent mixing of the two confinement 
effects. Nevertheless, eqs |5^, and ^ should be use- 
ful for the analysis of simulation data where real space 
distances are directly accessible. 

Experimentally, the microscopic deformationSpeAB|Only 
be measured via small angle neutron scattering.Eilc3 Un- 
fortunately, there seems to be no way to condense the 
structure functions eq |5^ which result from eqs |5^, |6^, 
and 1?^ for different strains into a single master plot. 
Figures ^ and || show a comparison for three character- 
istic values of A. Qualitatively, the results for the three 
models are quite similar. In particular, they do not pre- 
dict Lozenge-like patterns for the two-dimensional strufe. 
ture functions as they were observed by Straube ^ al.EHl 
In particular, we agree with Read and MrLeisha that 
the interpretation of Straube et al.E]E3c3 is based on 
an ad-hoc approximation in the calculation of structure 
functions from Model B. In principle, their alternative 
idea, to investigated the influence of dangling ends aa. 
the structure function within Model A and Model B,l1I 
can be easily extended to Model C. Judging from the 
small differences between the models (Figures ^ and ||) 
and the results in Ref.cfl, this would probably allow to ob- 
tain an excellent fit of the data and to correetJy account 
for the deformation dependence of the tubeE3 However, 
since the lozenge patterns were also observed in tri-block 
systems-^here only the central part of the chains was 
labeledpj dangling ends seem to be too simple an expla- 
nation. At present it is therefore unclear, if the lozenge 
patterns are a generic effect or iLliey are due to. 
artifacts such as chain scission.L£lEZl Simulationalil''^^' 
might help to clarify this point. 



IV. CONCLUSION 

In this paper we have presented theoretical consider- 
ations related to the entanglement problem in rubber- 
elastic polymer networks. More specifically, we have 
dealt with constrained fluctuation models in general and 
tube models in particular. The basic idea goes back to 
Edwards,El who argued that on a mean-field level dif- 
ferent parts of the network behave, as if they were em- 
bedded in a deformation-dependent elastic matrix which 



exerts restoring forces towards some rest positions. In 
the first part of our paper, we were able to show that the 
generalized Rouse modes of the corresponding phantom 
network without entanglement remain eigenmodes in the 
presence of the elastic matrix. In fact, the derivation of 
eq ^ which is the Hamiltonian of|-the exactly solvable 
Constrained Mode Model (CMM)JiJ provides a direct 
link between two diverging developments in the theory 
of polymer networks: the ideas of Edwards, Flory and 
others on the suppression of fluctuations due to entangle^. 
mentSp^d the considerations of Eichingerjid Graessley,EJ 
Markjla and others on the dynamics of (micro) phantom 
networks. An almost trivial conclusion from our theory 
is the observation, that it is not possible to estimate the 
entanglement effects from the knowledge of the absolute 
strength of the confining potentials alone. Required is 
rather the knowledge of the relative strength which in 
turn is a function of the network connectivity eq 

Unfortunately, it is difficult to exploit our formally ex- 
act solution of the constrained fluctuation model for arbi- 
trary connectivity, since it requires the eigenvalue spec- 
trum of the Kirchhoff matrix for randomly crosslinked 
networks. In the second part of the paper we have there- 
fore reexamined the idea of Heinrich and Straubecj to 
introcLuce entanglement effects into the Warner-Edwards 
modellj for linear, random paths through a polymer net- 
work, whose localization in space is modeled by a har- 
monic tube-like potential. In agreement with Heinrich. 
and Straube,E3 and with Rubinstein and Panyukovc2l 
we have argued that in contrast to confinement due to 
crosslinking, confinement due to entanglements is defor- 
mation dependent. Our treatment of the tube model dif- 
fers from previous attempts in that we explicitly consider 
the simultaneous presence of two different confining po- 
tentials. The effects are shown to be non-additive. From 
the solution of the generalized tube model we have ob- 
tained expressions for the microscopic deformations and 
macroscopic elastic properties which can be compared to 
experiments and simulations. 

While we believe to have made some progress, we do 
not claim to have solved the entanglement problem itself. 
For example, it remains to be shown how the geometrical 
tube constraint arises as a consequence of the topological 
constraints on the polymer conformations. But even on 
the level of the tube model, we are guilty of (at least) two 
possibly important omissions: (i) we have neglected fluc- 
tuations in the local strength of the confining potential 
and (ii) we have suppressed the anisotropic character of 
the chain motion parallel and perpendicular to the tube. 
In the absence of more elaborate tiieorieaj-computer sim- 
ulations along the lines of Refs.llJi23~E3'cII may present 
the best approach to a quantification of the importance 
of these effects. 
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FIG. 1. Excitation of constrained modes parallel and per- 
pendicular to the elongation at A = 1.5 as a function of the 
mode degree of confinement < 7 < 1. The dashed (dotted) 
lines show the predictions eq ^ of Model A and eq|| of Model 
B respectively for generalized Rouse modes of a phantom net- 
work with identical connectivity. The symbols represent the 
result of computer simulations of defect-free model polymer 
networks O The investigated modes are single-chain Rouse 
modes for network strands of length A'" « 1.25Ne. 
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FIG. 3. Plot of the parameters 2Ci and 2C2 of the 
Mooney-Rivlin equation / (A"^) = 2Ci + 2C2A"^ for the 
Rubinstein-Panyukov model (dotted) and the "double tube" 
model (solid). 
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FIG. 2. Langley plot of the shear modulus. The solid line 
corresponds to the "double tube" model, the dotted line to 
the Heinrich-Straube/ Rubinstein-Panyukov model and the 
dashed line to the phantom model. Ne represents the entan- 
glement length and Nc the crosslink length. 
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FIG. 4. Mooney-Rivlin representation of the reduced force 
for different values of $ (from top to bottom: the Phan- 
tom model (dashed line, "1? = 0), the "double tube" model 
(solid lines, $ — |,|,f) and the Heinrich-Straube/ Rubin- 
stein-Panyukov model (dotted line, "1> = 1)). 
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FIG. 5. Tube diameter dTc('l', A) = ((1 - $) + ^) ^ in 
parallel (upper curves) and perpendicular stretching direction 
for different elongation ratios A whereas $ can be expressed 
by the entanglement length Ne and the crosslink length Nc by 

,4 .2 

— — T usmg 



1+ 



= jf- . The dashed curve cor- 



responds to the Warner-Edwards model, i.e. drci^ ~ 0, A), 
the dotted curve corresponds to the Heinrich-Straube/ Ru- 
binstein-Panyukov model, i.e. drci^ ~ liA), and the solid 
line represents the "double tube" model with $ = |. 
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FIG. 6. Comparison of the universal scaling functions of 
eqs |5^, ^ for the Warner-Edwards model (dashed) and 
the Heinrich-Straube/ Rubinstein-Panyukov model (dotted) 

with y = . 
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FIG. 7. Kratky plots of the different structure factors in 
parallel and perpendicular stretching direction with ^ = 6: 
Warner-Edwards model (dashed), Heinrich-Straube/ Rubin- 
stein-Panyukov model (dotted line), "double tube" model 
with $ = I solid line). The upper curves correspond to the 
perpendicular stretching direction. 



14 



